Microscopic Analysis of Low-Frequency Flux Noise in YBa 2 Cu 3 07 Direct Current 
Superconducting Quantum Interference Devices 
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We use low-temperature scanning electron microscopy combined with SQUID detection of mag- 
netic flux to image vortices and to investigate low- frequency flux noise in YBa2Cu307 thin film 
SQUIDs. The low-frequency flux noise shows a nonlinear increase with magnetic cooling field up to 
60 fxT. This effect is explained by the surface potential barrier at the SQUID hole. By correlating 
flux noise data with the spatial distribution of vortices, we obtain information on spatial fluctuations 
of vortices on a microscopic scale, e.g. an average vortex hopping length of approximately 10 nm. 



oo 

(N 
O 

o 

Or 
C3 



-a 
c 

o 
o 



> 
m 

o 
o 

-I— > 
c3 



i 

C 

O 

o 



X 



PACS numbers: 68.37.Hk, 74.40.+k, 74.60.Ge, 74.76.Bz, 85.25.Dq 



The sensitivity of superconducting quantum interfer- 
ence devices (SQUIDs) based on high transition temper- 
ature T c superconductors is in most cases limited by low- 
frequency / flux noise, with the spectral density of mag- 
netic flux noise typically scaling as 5$ oc 1// (HEI- 
The main source of this excess noise is the thermally ac- 
tivated vortex motion in the high-T c thin films forming 
the body of the SQUID or the superconducting input cir- 
cuit 0. Low- frequency fluctuations induced by hopping 
of vortices are strongly dependent on the nature and spa- 
tial distribution of defects which act as pinning sites for 
vortices. Therefore, integral measurements of flux noise 
and electric transport properties give only very limited 
information on the underlying vortex dynamics. Hence, 
techniques offering spatial information are very useful to 
obtain a better understanding of the basic mechanisms 
of low-frequency flux noise generation. 

Over the last two decades, low-temperature scanning 
electron microscopy (LTSEM) has proven to be very suc- 
cessful in getting local information on properties of su- 
perconducting thin films and Josephson junctions, such 
as the spatial distribution of critical temperature T c or 
critical current density j c and of Josephson vortices in 
long junctions 0. Recently, this method was successfully 
extended to the direct imaging of vortices in YBa2Cu307 
(YBCO) SQUIDs 4]. Together with a simultaneous mea- 
surement of the low-frequency flux noise of the SQUIDs 
one obtains both images of the spatial distribution of 
vortices and local information on noise properties of the 
devices under investigation^- 

In this paper we investigate YBCO dc SQUIDs by LT- 
SEM to study the increase of magnetic flux noise with 
magnetic cooling field and to analyze the mechanism 
of low-frequency flux noise generation on a microscopic 
scale. 

In our experiments, we use epitaxially grown YBCO 
thin film dc SQUIDs of the so-called washer design [see 
Fig. Ufa)] • The c-axis oriented YBCO film has a thickness 
of 80 nm; the washer size is 120 x 305 /im 2 , with a 100 
fj.ui long and 4 /jm wide slit. The 1 ^m wide Josephson 
junctions are formed by a 24° symmetric grain boundary 
in the underlying SrTi0 3 substrate. The YBCO SQUIDs 
are mounted on a magnetically shielded, liquid nitrogen 
cooled cryostage of a SEM yfl and read out by a stan- 
dard flux-locked loop (FLL) with 3.125 kHz bias current 
reversal to eliminate 1// noise due to fluctuations in the 
critical current I c of the Josephson junctions. 
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FIG. 1: (a) SQUID washer design; dotted line indicates grain 
boundary, (b) Coupling strength cf> r vs. radial distance r from 
the slit [numerical calculation for design in (a)]. Inset: Mag- 
netic flux $/"!>o(r) coupled by a single vortex into the SQUID 
hole; numerical calculation for design in (a) [squares] and for 
a circular washer [disks]; the line is obtained from analytical 
calculation 6] for a circular washer, (c) 5$ images of the 
SQUID washer for two different cooling fields Bq (T"=77K). 



For the spatially resolved measurements, the electron 
beam is used as a local perturbation which induces an 
increase in temperature 8T(x — Xq, y — yo) on the sample 
surface (in the x-y-plane) centered around the beam spot 
position (xo,yo). The length scale for the spatial decay 
of the thermal perturbation is set by the beam-electron 
range R w 0.5 /xm for a typical beam voltage Vt, = 10 
kVj3(- This gives a maximum increase in beam induced 
temperature AT w 1 K at (xo,yo) for a typical beam 
current lb — 7 nA. So-called 5&(x , y ) images are ob- 
tained by recording the e-beam induced flux change S$ 
in the SQUID as a function of the e-beam coordinates 
(xo,yo). To improve the signal to noise ratio, we use a 
beam-blanking unit operating at typically 5 kHz and the 
output signal of the FLL, i.e., the e-beam induced flux 
change in the SQUID, is lock-in detected. Additionally, 
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the time trace or the power spectrum of the FLL output 
signal can be recorded by a signal analyzer. 

The mechanism of imaging of vortices is explained in 
Ref. H and can be briefly described as follows: The e- 
beam-induced local increase in temperature induces a lo- 
cal increase in the London penetration depth A^. Hence, 
the screening currents circulating around a vortex are 
spatially extended due to e-beam irradiation. If the e- 
beam is scanned across a vortex this vortex is virtually 
dragged along with the beam, i.e., it is displaced by some 
distance 5r, if the beam spot is within the radial dis- 
tance R from the vortex. This displacement changes the 
amount of stray magnetic flux that a vortex couples to 
the SQUID hole. Hence, scanning across a vortex induces 
a negative (positive) flux change J$ in the SQUID if the 
vortex is moved away from (towards) the SQUID hole. 
Figure ^c) shows two examples of <5$ images with vor- 
tices appearing as pairs of positive (bright) and negative 
(dark) signals. These images were taken after cooling 
from above T c to T = 77 K in a static magnetic field Bo 
— 60 /iT and 5 /iT. The maximum displacement Ar of a 
vortex is of the order of the beam-induced change in Al 
(typically w 20 nm in our experiments at 77 K) 5j. The 
displacement Ar induces a signal <5$ = ^(r)x Ar, with 
r being the radial distance of the vortex from the SQUID 
hole. 

For a given washer geometry, the coupling strength 
(j)r(r) = is a function of the vortex position only. 

Therefore, it can only be obtained from a spatially re- 
solved measurement. <j) r plays an important role for the 
signal noise as well: The flux noise S$i from a moving 
vortex is 5$i = </>r( r «) x ^rit where S r i is the spectral 
density of radial motion of the vortex at distance r-j p| . 

By numerical simulation^ one can calculate the 
amount of stray flux $ that a single vortex couples into 
the SQUID hole. For a given SQUID design, this value 
to a good approximation depends for vortices to the left 
or right of the SQUID hole [see Fig^a)] omv on the 
distance r from the SQUID hole. The dependence <&(r) 
(normalized to the flux quantum $o) is shown in the inset 
of Fig.^b). We note that analytical calculations of $(r) 
for circular washers U are in excellent agreement with our 
numerical simulation results [c.f. inset of Fig. Gib)]. We 
find that $(r) strongly depends on the washer geometry; 
e.g. $(r) falls off with increasing r much more rapidly for 
a circular washer as compared to a square washer. Gener- 
ally, <I>(r) quickly falls off from $o( r = 0) with increasing 
r and approaches zero at the outer washer edge. The 
derivative of this function, the coupling strength <?V(r), 
shows a similar behavior [see Fig. ^b)] . This has two 
important effects on vortices located close to the SQUID 
hole: On the one hand, the signal contrast in the image 
is high, since <5i> oc <j) r . On the other hand, the noise sig- 
nal produced by thermally activated motion of the vortex 
contributes strongly to the total noise, since S<& oc (f) r . 

Vortices in superconducting thin films show thermally 
activated hopping between two (or more) pinning sites. 
Due to the small hopping length, the hopping process 
cannot be resolved directly by LTSEM, but it appears as 
a random telegraph signal (RTS) in the SQUID response 
due to a change in the magnetic flux coupled to the 
SQUID hole. Sometimes, the SQUID signal is dominated 
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FIG. 2: Rms flux noise of YBCO dc SQUID at 77 K showing 
Lorentzian noise due to a single dominating fluctuator (a) and 
low- frequency l//-noise (b). Inset: time trace of the SQUID 
output with a single dominating fluctuator. 



by a single fluctuator. In this case the RTS is clearly 
visible [see inset in Fig. [2], producing a Lorentzian- 
shaped noise spectrum^ [see trace (a) in Fig. 0- An 
uncorrelated superposition of many RTSs leads to a 1//- 
shaped noise spectrum, provided the energies of the po- 
tential barriers between the pinning sites are distributed 
uniformly |lflj [see trace (b) in Fig. E). 

For uncorrelated fluctuations of N vortices trapped in 
the washer, the total noise is 

N 
i=l 

where S r is the average spectral displacement noise power 
and K = X)i^r( r i) i s the sum of the squared coupling 
strengths. With the vortex positions rj, obtained from 
the (5<I> images (for different Bq, i.e. different N) and 
the calculated function (f) r (r) we determine K(N) [see 
Fig. Ej. Surprisingly, we find K oc N 1 - 5 instead of 
K hom oc N, which is expected for a homogeneous distri- 
bution of vortices, and which we calculate by replacing 
the sum with the integral N x J 4>1{r)dr. Furthermore, 
we find K hom > K for all N [see dashed line in Fig. El- 

The reason for this deviation is a superproportional 
increase of the density of vortices close to the SQUID 
hole with increasing Bq, as imaged by LTSEM [see inset 
of Fig. |3] ■ This result can be explained by the forma- 
tion of a surface barrier for vortex exit at the washer 
edges: For low magnetic fields, i.e. few vortices in 
the washer structure, the vortices are far apart from 
each other. In this case the repulsive surface barrier 
potential [Tlj is large compared to the repulsive vortex- 
vortex interaction 1 1 2|. preventing formation of vortices 
close to the SQUID hole during the cooling process. For 
higher fields, i.e. increased vortex density, the vortex- 
vortex repulsion becomes stronger, thus pushing vortices 
closer to the SQUID hole. 

With S$ oc K and the scaling K oc iV 15 shown in 
Fig. |3 one expects S$ to increase stronger than linear 
with the number of vortices N or magnetic field. Figure 
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FIG. 3: Sum of squared coupling strengths K vs. number 
of vortices N . Dashed line: homogeneous vortex distribution. 
Inset: fraction of vortices found at a distance < 6 /mi from 
the slit edge for different cooling fields Bo- 
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FIG. 4: Spectral density of flux noise S$(lHz) vs. number 
of vortices N. Dashed line: homogeneous vortex distribu- 
tion. Inset: Average spectral density of radial motion SV(lHz) 
vs. N, calculated from K and S$(lHz). 



Elshows data for S$(lHz) vs. N. Although S$(iV) scat- 
ters quite strongly, it can be clearly seen that it increases 
more like oc TV 1,5 than oc N, in contrast to previous mea- 
surements on larger SQUID washers [J. By measuring 



5$(lHz) and determining K from the 5$ pictures, one 
can calculate S r (l Hz) using S<& = S r X K [see inset of 
Fig.|H. Hz) is found to be 0.2-1.5 nm 2 /Hz. For 

a homogeneous vortex distribution, assuming S r = 0.7 
nm 2 /Hz, one would get Sl° m oc N with Sl° m > 5$ in 
the investigated field range [see dashed line in Fig. 2] ■ 

From these data it is po ssible to estimate the average 
vortex hopping length Ax using a model based on the 
following assumptions: (i) All vortices are situated in 
a symmetric double well potential with fixed distance 
Ax r of local minima and variable barrier energy E. (ii) 
The barrier energies E are distributed uniformly, such 
that the highest (lowest) characteristic frequency f„[To| 
is 1000 Hz (0.01 Hz). This condition provides an 1//- 
shaped spectrum. Now Ax r can be calculated using the 
flux noise theory developed in Refs. [T0| and [J. The 
result is Ax r « (3-8) nm. This calculation gives the 
hopping length in radial direction. Considering hopping 
in two dimensions, these numbers have to be multiplied 
by the factor ir/2 to obtain the overall average hopping 
length Ax « (5 - 13) nm. Although this model uses 
some rather simple approximations it should produce the 
correct order of magnitude of the vortex hopping length. 

In conclusion, we combined vortex imaging by LTSEM 
with measurements of low-frequency flux noise in YBCO 
dc SQUIDs to obtain information on spatial fluctuations 
of vortices on a microscopic scale. We have shown that a 
nonlinear increase of the flux noise with magnetic cooling 
field can be explained by a surface barrier effect. Our 
results clearly show that either strong pinning or avoding 
of vortices trapped close to the SQUID hole is essential 
for reducing low- frequency flux noise. Furthermore, our 
analysis yields an average vortex hopping length around 
10 nm for our films. The nature of defects producing 
pinning sites with hopping lengths on this scale still has 
to be revealed. 
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